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We study the dependence on the spatial dimensionality of different quantities relevant in the description of the 
Anderson transition by combining numerical calculations in a 3 < d < 6 disordered tight binding model with 
theoretical arguments. Our results indicate that, in agreement with the one parameter scaling theory, the upper 
critical dimension for localization is infinity. Typical properties of the spectral correlations at the Anderson 
transition such as level repulsion or a linear number variance are still present in higher dimensions though 
eigenvalues correlations get weaker as the dimensionality of the space increases. It is argued that such a critical 
behavior can be traced back to the exponential decay of the two-level correlation function in a certain range of 
eigenvalue separations. We also discuss to what extent different effective random matrix models proposed in 
the literature to describe the Anderson transition provide an accurate picture of this phenomenon. Finally, we 
study the effect of a random flux in our results. 
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I. INTRODUCTION 

After almost fifty years of the landmark paper by 
Anderson 1 about localization, the study of the properties of 
a quantum particle in a random potential is still one of the 
central problems of modern condensed matter physics. 

In the the early days of localization theory, research was 
largely focused on the determination of the critical disorder at 
which the the metal-insulator transition -also referred as the 
Anderson transition- occurs as a function of the connectiv- 
ity of the lattice. In the original Anderson's paper this was 
achieved by looking at the limits of applicability of a loca- 
tor expansion^ Later on, a more refined estimation based 
on the solution of a self-consistent equation 3 provided with a 
similar answer. The self-consistent method is only exact in 
the case of Cayley tree but it is believed to be accurate if the 
spatial dimensionality is large enough. We note that in the lo- 
cator expansion 1 the metal-insulator transition is induced by 
increasing the hopping amplitude of an initially localized par- 
ticle. 

In the seventies the application of ideas and techniques of 
the theory of phase transitions such as scaling and the renor- 
malization group 4 opened new ways to tackle the localiza- 
tion problem especially in low dimensional systems. These 
progress led eventually to the proposal of the one parameter 
scaling theory 5 which, despite being still under debate, has be- 
come the 'standard' theory of localization. In the one parame- 
ter scaling theory, localization in a given disordered sample is 
described by the dimensionless conductance g. This quantity 6 
is defined either as the sensitivity of a given quantum spec- 
trum to a change of boundary conditions in units of the mean 
level spacing A or as g = E c /A where E c , the Thouless en- 
ergy, is an energy scale related to the classical diffusion time 
to cross the sample. The dimensionless conductance g is sen- 
sitive to localization effects. In a metal (insulator), it increases 
(decreases) monotonically with the system size, L. 



Under the assumption that the dimensionless conductance 
is only a function of the system size and by using simple 
scaling arguments, the one parameter scaling theory predicts 
that the metal-insulator transition is characterized by a scale 
invariant dimensionless conductance g = g c . The lowest 
dimension in which the metal-insulator transition occurs is 
d > 2. In two and lower dimensions destructive interference 
caused by backscattering produces exponential localization of 
the eigenstates in real space for any amount of disorder in the 
limit L — > oo. In this picture, the Anderson transition is con- 
sidered as a standard second order transition with critical ex- 
ponents s , v that control how the conductivity a oc | W — W c | s 
vanishes or the localization length £ oc \W C — diverges 
as the critical disorder W c is approached. 

In d = 2 + e (e <C 1) the transition occurs in the weak 
disorder region and consequently an analytical treatment is 
possible. Diagrammatic perturbation theory and field theory 
techniques 4,7 predict that v ~ 1/e and W c oc e. By con- 
trast, the critical exponent associated to the Cayley tree, which 
should be close to that of a disordered conductor in d 3> 2 di- 
mensions, is v = 1/2. 8 In the context of second order phase 
transitions this value corresponds with the upper critical di- 
mension d u , namely, for d > d u fluctuations are irrelevant 
and the mean field approximation become exact. For the lo- 
calization problem different values d u = 4, 6, 8, oo of the up- 
per critical dimension have been reported. 9 The results of this 
paper discard d u — 4, 6 and indicate that d u — » oo is the up- 
per critical dimension. However we would like to point that 
the exact significance of the upper critical dimension for lo- 
calization is unclear. It is not known what fluctuations are 
suppressed at the upper critical dimension and to what extent 
spectral or transport properties at criticality are affected. 

The Anderson transition in a disordered conductor is a con- 
sequence of a highly non trivial interplay between quantum 
destructive interference effects and quantum tunneling. In 
low dimensions, d ~ 2, weak quantum destructive interfer- 
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ence effects induce the Anderson transition. Analytical results 
are available based on perturbation theory around the metal- 
lic state. 4 ' 7 In high dimensions, d 3> 2, quantum tunneling is 
dominant and the locator expansion 1 or the the self-consistent 
formalism 3 can be utilized to describe the transition. We note 
that in these papers corrections due to interference of different 
paths are neglected. 

The progress in numerical calculations during the 
last twenty years has increased dramatically our 
knowledge 10 i n ' 12 - 13 i 14 of the metal-insulator transition 
specially in intermediate dimensions such as d = 3,4 for 
which a rigorous analytical treatment is not available. Below 
we cite a few of its most relevant results. 

It was verified that for a disorder strength below the crit- 
ical one, the system has a mobility edge at a certain energy 
which separates localized from delocalized states. 10 Its posi- 
tion moves away from the band center as the disorder is de- 
creased. Delocalized eigenstates, typical of a metal, are ex- 
tended through the sample and the level statistics agree with 
the random matrix prediction 15 for the appropriate symme- 
try. The spectral correlations at the Anderson transition, usu- 
ally referred to as critical statistics 1 1,ie , are scale invariant and 
intermediate between the prediction for a metal and for an 
insulator. 11 By scale invariant we mean any spectral correla- 
tor utilized to describe the spectral properties of the disordered 
Hamiltonian does not depend on the system size. 

Eigenfunctions at the Anderson transition are 
multifractals 1013 (for a review se o 17 ' 18 ), namely, their mo- 
ments present an anomalous scaling, V q = J d d r\ijj(r)\ 2q oc 
L-Dqiq- 1 ) w ith respect to the sample size L, where D q is a 
set of exponents describing the Anderson transition. 

The main features of the Anderson transition only de- 
pend on the dimensionality of the space and the universality 
clas a 20 i 21 , namely, the presence or not of a magnetic field (or 
other time reversal breaking mechanism) or a spin-orbit inter- 
action. The dependence with the universality class diminishes 
as the spatial dimension increases. It has also been reported 
that certain spectral correlators at the Anderson transition are 
sensitive to different boundary conditions* 2 ^ 

All of these numerical findings are compatible with the one 
parameter scaling theory. The applicability of the e-expansion 
(d = 2 + e) is by contrast much more restricted. A naive ex- 
trapolation to e = 1, 2 yields vj,o ~ 1/e = 1, v±r> — 1/2 thus 
suggesting that the upper critical dimension is four. However 
numerical calculations 12 ' 23 show undoubtedly that v?,d ^3/2 
and v^ei ~ 1. Similarly, up to d = 4, the self-consistent the- 
ory overestimate by more than a factor two the value of the 
critical disorder at which the Anderson transition occurs. This 
suggests none of the theories traditionally utilized to describe 
the metal-insulator transition can be really extrapolated to the 
physically relevant case of d — 3. In order to make progress in 
this difficult problem a new basis for the study of the Ander- 
son transition in any dimension is necessary. In this paper we 
have a more modest goal: a detailed exploration of the depen- 
dence of different quantities defining the Anderson transition 
on the spatial dimensionality. 

We propose simple relations that describe how the param- 
eters defining the Anderson transition depend on the dimen- 



sionality of the space. It is argued that the upper critical di- 
mension must be infinite. Our results are supported by numer- 
ical evidence from a disordered Anderson model in a hyper- 
cubic lattice in 3 < d < 6. This is the first time than the 
Anderson transition in d = 5, 6 is investigated numerically in 
the literature (for some recent results in a small asymmetric 
lattice in d = 5 we refer to Refill). 

The organization of the paper is as follows. In section II, we 
introduce the model to be studied, explain the technical details 
of the numerical simulation, locate the mobility edge in differ- 
ent dimensions and investigate how the critical exponent and 
the critical disorder depend on the spatial dimensionality. In 
Sec. EH, spectral correlation at the Anderson transition are in- 
vestigated, three different regions of the two level correlation 
function are distinguished, we also study the dependence of 
the slope of the number variance and the asymptotic decay of 
the level spacing distribution with the spatial dimensionality. 
It is also discussed the range of validity of certain phenomeno- 
logical models commonly used in the literature to describe the 
spectral correlations at the Anderson transition. Finally, we 
examine the effect of a magnetic flux in our results. 



II. CRITICAL DISORDER, CRITICAL EXPONENTS AND 
UPPER CRITICAL DIMENSION 

In this section we determine the critical disorder and criti- 
cal exponents for different dimensions and then discuss their 
dependence with the spatial dimensionality. 



A. The model: Technical details 

Our starting point is the standard tight-binding Anderson 
model on a hyper-cubic L d lattice with d = 3, . . . , 6 

U = ^ e t a i a i + X! t ij a i a 3 j W 
i ij 

where the operator ai{a\) destroys (creates) an electron at the 
ith site of the lattice and tij is the hopping integral between 
sites i and j which is non zero only for nearest neighbors. In 
the following we take t ij = 1 for i , j and the lattice constant 
equal to unity, which sets the energy and length units, respec- 
tively. The uncorrelated random energies ej are distributed 
with constant probability within the interval (— W/2, W/2), 
where W denotes the strength of the disorder and hard-wall 
boundary conditions are imposed in all directions. 

In order to proceed we compute eigenvalues of the Hamil- 
tonian Eq. (Q3 for different volumes and disorders by using 
techniques for large sparse matrices, in particular a Lanczos 
tridiagonalizaton without reorthogonalization method^ 5 - We 
restrict ourselves to a small energy window (—2,2) around 
the center of the band. Calculations have been carried out in 
samples of sizes up to L = 30 for d = 3, 12 (d = 4), 10 
(d = 5) and 7 (d — 6). The number of random realizations is 
such that for a given triad of {d, L, W} the number of eigen- 
values obtained is at least 3 x 10 5 . In order to study the level 
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FIG. 1 : Scaling variable rj as a function of disorder W for different 
volumes in d = 5 (left panel) and d — 6 (right panel). Energies 
in the interval ( — 2, 2) have been considered. The system undergoes 
an Anderson transition at W c = 51.4 and 74.5 for d — 5 and 6, 
respectively (see text). 



statistics around the mobility edge more accurately, this num- 
ber was increased to 20 x 10 6 at the critical disorder. 

Eigenvalues thus obtained are appropriately unfolded, i.e, 
they were rescaled so that the spectral density on a spectral 
window comprising several level spacings is unity. 



B. Location of W c (d) and v{d) 



sition is found in such a high dimensional disordered system. 
For a precise determination of the critical disorder W c and the 
critical exponent v we look at the correlation length near W c 



t(W) = So\W-W c 



(3) 



where £o is a constant. The numerical values of W c and v 
are obtained by expressing r)(L, W) = f[L/£(W)] and then 
performing an expansion around the critical point 



rj(L, W)^r lc + J2c n (W - W c ) n V 



ijv 



(4) 



In practice, we have truncated the series at n = 4. For each di- 
mension (d — 5, 6) we have performed a statistical analysis of 
the data in the windows shown in Fig. 1 with the Levenberg- 
Marquardt method for nonlinear least-squares models. The 
most likely fit is determined by minimizing the \ 2 statistics 
of the fitting function Q. We found the following critical dis- 
orders W c = 51.4 ± 0.4 in d = 5 and W c = 74.5 ± 0.7 
in d = 6, and the corresponding critical exponents are equal 
to v = 0.84 ± 0.06 and v = 0.78 ± 0.06, respectively. A 
similar analysis for the d = 3 and d = 4 systems results in 
W c = 15.22 ± 0.08 and v = 1.52 ± 0.06 for d = 3, and 
W c = 29.8 ± 0.2 and v = 1.03 ± 0.07 for d = 4. We note 
that in the d = 3 case, the deviation of W c from the accepted 
value W c ~ 16.5 is due to the utilization of rigid boundary 
conditions. See Fig. 2 for a plot of W c and v as a function of 
the spatial dimensionality. 



Our first task is to find out the critical disorder W c and the 
critical exponent v for different dimensions in the small spec- 
tral window (—2, 2) around to the origin. In order to proceed 
we determine the location of the mobility edge close to the 
band center by using the finite size scaling method. 1 1 First we 
evaluate a certain spectral correlator for different sizes L and 
disorder strengths W. Then we locate the mobility edge by 
finding the disorder W c such that the spectral correlator ana- 
lyzed becomes size-independent. In our case we investigate 
the level spacing distribution P(s) (probability of finding two 
neighboring eigenvalues at a distance s = (ej+i — Ci) /A, with 
A being the local mean level spacing). The scaling behav- 
ior of P(s) is examined through the following function of its 
variance^ 

r](L, W) = [var(s) - var WD ]/[ var P - var W o] , (2) 

which describes the relative deviation of var(s) from the 
Wigner-Dyson (WD) limit. In Eq. © var(s) = (s 2 ) - (s) 2 , 
where (. . .) denotes spectral and ensemble averaging, and 
varwD = 0.286 and varp = 1 are the variances of WD and 
Poisson statistics, respectively. Hence rj = 1(0) for an insula- 
tor (metal). Any other intermediate value of r] in the thermo- 
dynamic limit is an indication of a mobility edge. 

In Fig. 1 we plot the W dependence of rj for different sys- 
tem sizes in d = 5 (left panel) and d = 6 (right panel). The 
critical disorder W = W c signaling the Anderson transition 
corresponds with the point for which r) is independent of L. 
For a weaker (stronger) disorder, r] tends to the metallic (insu- 
lator) prediction. This is the first time that an Anderson tran- 



C. Theoretical analysis of W c (d) and v(d) 

In certain limiting cases W c and v are known analytically. 
For instance, if effects related to interference among differ- 
ent paths are neglected 3 , the standard tight binding Anderson 
model is effectively defined on a Cayley tree and v = 1/2. On 
the other hand if only interference corrections to the metallic 
limit are considered then v = l/(d — 2). 7 The former predic- 
tion is supposed to be approximately valid for d 3> 2 and the 
latter for d = 2 + e and e <C 1. From the above numerical 
results it is clear that none of these limits is appropriate in the 
range of intermediate dimensions of interest. Additionally, it 
is believed 3 that corrections to the v = 1/2 result should go as 
~ 1/d since this is the dependence on the spatial dimension- 
ality of the neglected diagrams describing interference effects. 
Combining these two facts we propose that 



for all dimensions. As is shown in Fig. 2 (right panel), this re- 
lation verifies all limiting cases and reproduce the numerical 
results accurately. According to the above relation, the up- 
per critical dimension for localization is infinite. This result 
is fully supported by the analysis of the spectral correlations 
(see next section). Moreover, in a recent paper 27 it has been 
proved rigorously that the level statistics of a disordered sys- 
tem in a Cayley tree (y = 1/2) is Poisson as for an insulator. 
As was mentioned previously, the Cayley tree represents to 
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FIG. 2: W c (left panel) and v (right panel) as a function of the di- 
mensionality d. The simple interpolating formula (solid line) Eq. (5) 
describes quite well the numerical results for the critical exponents 
v. 

a d-dimensional conductor where all interference effects be- 
tween different paths are neglected. It is thus supposed to be 
an accurate description of a disordered conductor only in the 
limit d > 2. 

On the other hand, if the one parameter scaling theory is 
valid, quantum diffusion never stops (see next section) for any 
finite dimension. Level repulsion typical of a metal will be 
present in any finite dimension so level statistics at criticality 
can be that of an insulator only in the d — > oo limit. But this 
precisely the result for the Cayley tree 27 which corresponds 
with the upper critical dimension for localization. It is thus 
clear that the upper critical dimension must be infinity. 

A similar analysis can be carried out for the critical dis- 
order W c . The original estimation of Anderson^, W c — 
4K \n(W c /2) (where K is the connective constant which is 
a little bit less than the number of nearest neighbors minus 
one) greatly overestimates W c . This is hardly surprising since 
the Anderson's calculation involves crude approximations and 
consequently should be considered as an order of magnitude 
estimation rather than an accurate prediction. For instance, 
deviations coming from interference effects are neglected in 
this scheme. Roughly speaking they tend to reduce W c by an 
amount of order W c /d. 

In the opposite limit, d = 2 + e, simple perturbation theory 7 
yields W c oc d — 2. The discrepancy observed with the an- 
alytical results in the limit of high dimensionality prevent us 
from proposing an interpolating relation as in the case of the 
critical exponent. However we have noticed that a much better 
agreement with the numerical results is achieved if an effec- 
tive connective constant K e s — K/2 is utilized (solid line 
in left panel of Fig. 2). Furthermore, the remaining devia- 
tion gets smaller as the spatial dimensionality increases thus 
suggesting that it may be produced by destructive interference 
effects (~ W c /d). 



III. LEVEL STATISTICS 

In this section we investigate the level statistics at the An- 
derson transition. We shall mainly focus on its dependence 



on spatial dimensionality and the exact functional form of the 
two-level correlation function (TLCF). 



A. Theoretical analysis of R2 (s) 

Our starting point is the connected TLCF, 

Ms) = j^^ e - w / 2 M e + w / 2 )> > < 6 > 

where p(e) is the density of states at energy e, ( } denotes 
averaging over disorder realizations and s = u>/ A where A = 
1/L d {p(e)) in the mean level spacing. Once the spectrum has 
been unfolded J?2(s) can be simply written as 

R 2 (s) = S(s)+J2p(n;s) , (7) 

n 

where p(n; s) is the distribution of distances s n between 
n other energy levels and 8{s) describes self-correlation of 
levels. 15 In numerical computations we use Eq. © since it 
gives much more accurate results than Eq. ©. 

According to the one parameter scaling theory, the spectral 
properties depend on the dimensionless conductance g which 
is a function of the system size L only. In a metal g — > 00 for 
L — > 00, the Hamiltonian can be accurately approximated by 
a random matrix with the appropriate symmetry and Wigner- 
Dyson statistics applies. 15 For instance, for broken time re- 
versal invariance, R^^s) = S(s) + 1 — sl " . In an insu- 
lator, eigenvalues uncorrelated, Poisson statistics applies and 
R 2 (s)=6(s). 

Right at the Anderson transition, the dimensionless con- 
ductance g = g c is size independent and level statistics are 
supposed to be universal and intermediate between Wigner- 
Dyson and Poisson statistics. Unfortunately there are few ana- 
lytical results for the TLCF at criticality. In the d > 2 region 4 , 
g c ~ — 2) ^> 1. The TLCF, can only evaluated ex- 
plicitly in the limit — » 2 28 and g ^> 1 where i?2(s) ~ s (time- 
reversal), for s <C g c , as for a metal, and i?2(s) ~ cr As / 9c for 
s 3> g c , with A being a factor of order unity. The Anderson 
transition is thus characterized by level repulsion combined 
with an exponential decay of the TLCF. 

In higher dimensions the exact form of the TLCF is not 
known. However, we note that the scale invariance of the 
spectral correlations at the Anderson transition restricts the 
decay of the TLCF in the s 3> g c region to be power-law or 
exponential 19 . Our numerical results (see Fig. 3) for d > 3 
support also this picture. 

The limit of long times and small energy differences s -c 
g c , s <§; g c is well understood in high dimensions as well. 
Level repulsion of neighboring eigenvalues i?2(s) °c s, typi- 
cal of a metal, should be a generic feature in any dimension. 
According to the one parameter scaling theory, the averaged 
moments of the particle position at the Anderson transition 
increase asymptotically t — > 00 as (r(t) 2m ) ~ ^m/d wnere 
to is a positive integer. As the spatial dimensionality d in- 
creases the diffusion is slowed down but it never stops even 
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for long times. This is an indication that the spectral correla- 
tions for sufficiently small energy intervals are similar to those 
of a metal and, as a consequence, i?2 (s) ~ s for s <§C g c (see 
Fig. 3). 

Finite size effects modify the TLCF in the critical region. 29 
In any finite system at criticality the localization length £ oc 
\E C — E\~ v (E c is the location in energies of the mobil- 
ity edge) is finite and the dimensionless conductance is not, 
strictly speaking, scale invariant, g(L^) = g c [l + (L^/L) 1 ^] 
where is the localization length for a given E ~ E c , 
As a consequence 29 , the TLCF develops a power-law tail, 
R$ il (s) oc s^ 2 with 7=1- l/(vd) for s > A e /A, where 
Aj is the mean level spacing in a localization volume £ d . This 
tail is not related with the properties of the critical point but 
rather with how the system approaches to it. In d = 2 + e, 
v = 1/e > 1 and Rf ll (s) ~ 1/s. 

As a summary, we can distinguish three different regions in 
the critical TLCF, for s <C g c , R2(s) oc s, for s 3> g c , R2(s) 
decays exponentially. For s > A^/A decays as power-law 
due to finite size effects. In order to observe the exponential 
decay related to the critical point our system size must be such 
that g c > A^/A. Finally we note that the exact dependence 
of g c on the spatial dimensionality it is not known. We are 
only aware of the prediction of Vollhardt and Wolfle 7 by using 
a self-consistent diagrammatic theory valid for 4 > d > 2, 
9c(d) = c d /(d - 2) with c d = (2/n)[S d /(2n) d } and S d the 
surface of a ri-dimensional sphere of radius unity. In principle 
it should be accurate only for d > 2 though it is unclear its 
exact range of validity. 



B. Numerical analysis of f?2(s) 

After the theoretical analysis we are now ready to present 
our numerical results for i?2(s) at the Anderson transition in 
d = 3 — 6 dimensions. Our motivation is to study the ex- 
istence and extension of the three regions introduced above: 
level repulsion, power-law and exponential decay. Indeed our 
numerical results clearly show these three regimes in all di- 
mensions d = 3 — 6 investigated. 

We have first verified (not shown) that for sufficiently large 
s, i?2(s) ~ 1/s 7 - The numerical value of the exponent 7 
was in full agreement with the theoretical prediction 7=1 — 

(yd)- 1 . 

Then we investigate to what extent level repulsion typical of 
the Anderson transition in d = 3, 4 is still present in higher di- 
mensions. As is observed in Fig. 3, left, for sufficiently small 
s, i?2(s) ~ s for all dimensions studied. The solid lines lines 
are linear fits of the form R2 (s) = C + Ds with fitting param- 
eters D = 6.6 ± 0.8 for d = 3, 15.0 ± 1.2 (d = 4), 101 ± 5 
(d — 5), and 373 ± 32 (d = 6). The parameter C is equal to 
zero within the error bars in all cases. This is consistent with 
the prediction of the one parameter scaling theory that quan- 
tum diffusion never stops. However the range in which level 
repulsion is observed decreases dramatically with the spatial 
dimensionality thus suggesting that the critical conductance 
g c also decreases rapidly with the dimension. It is hard to give 
a more quantitative prediction of g c as a function of the spatial 
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FIG. 3: i?2(s) for d = 3 — 6 at the Anderson transition for different 
s ranges. In the left panel we look at the region of small s where 
level repulsion is still observed. As the dimensionality increases the 
Anderson transition occurs for stronger disorder and the region of 
level repulsion is smaller. In the right panel it is shown the window 
of a in which exponential decay is observed. Such decay is responsi- 
ble of typical features of the Anderson transition as a linear number 
variance or a scale invariance spectrum. For the sake of clarity we 
have removed the power-law contribution i?2(s) oc 1/s 7 . It is well 
established that this term does not really describe the properties at 
the Anderson transition but rather how the system approaches to it. 
Moreover its contribution to the number variance and other spectral 
correlators is negligible with respect to the exponential contribution. 



dimensionality: the estimation of Vollhardt and Wolfle 7 men- 
tioned previously fails for d > 3. Another option is to extrap- 
olate the result in the diffusive regime 30 i?2(s) ~ s(l+a d /g 2 ) 
for s g c to the critical one. However the geometrical coef- 
ficient a d diverges also for d > 3. 

Our numerical results (see Fig. 3, right) show that for larger 
spectral separations s > g c , the linear repulsion is replaced by 
an exponential decay. The solid lines correspond to a linear 
fit ln[l — i?2(s)] = C — Ds with fitting parameters D = 
3.7 ± 0.1 for d = 3, 4.7 ± 0.1 (d = 4), 5.6 ± 0.5 (d = 5), 
and 9.5 ± 0.2 {d — 6). The maximum value of s plotted 
was chosen attending to technical criteria. For larger values 
of s, 1 — i?2(s) fluctuates around zero thus suggesting that the 
maximum precision of the computer has been reached. 

We note that such exponential decay has already observed 
in certain one-dimensional disordered systems with long 
range hopping 13 and in phenomenological short-range plasma 
models 34 whose spectral properties are strikingly similar to 
those of a disordered system with short range hopping at the 
Anderson transition. In these one-dimensional systems it can 
be proved analytically that R c (s) ~ e~ As 1 9 where A is a con- 
stant of order unity. It is thus tempting to speculate that in our 
case g c ~ 1/D with D the fitting parameter above. However 
a deeper analytical knowledge about the Anderson transition 
is needed to discard that additional geometrical factors (as a d 
above) enter in the exponent of i?2(s) making thus less evi- 
dent the relation between g c and D. 
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FIG. 4: Number variance E 2 (i) (left panel) and P(s) (right panel) 
at the Anderson transition for d — 3 — 6. The system tends to the 
Poisson limit (P) as the dimension is increased. WD denotes the 
Wigner-Dyson distribution. 



C. Spectral correlators 

Level statistics at the Anderson transition are usually in- 
vestigated by computing certain spectral correlators from the 
TLCF or higher rt-level correlation functions. The level spac- 
ing distribution P(s) is a popular choice to study the corre- 
lations of eigenvalues separated short distances of order the 
mean level spacing. On the other hand, the number variance 
Y?(£) = ((Ni - (N e )) 2 ) (N e is the number of eigenvalues 
in an interval of length £) provides useful information about 
spectral correlations for distances much larger than the mean 
level spacing. 

Numerical calculations in d < 4 at the Anderson transition 
have found that P(s) — > for s — > 0, as in a metal. However 
the number variance is asymptotically linear E 2 (I) ~ x&> as 
in an insulator but with a slope x < 1- The origin of this linear 
behavior can be explained heuristically 31 by using the one pa- 
rameter scaling theory and making the plausible approxima- 
tion that eigenvalues interact only if their separation (in units 
of the mean level spacing) is smaller than g c . In the critical re- 
gion it is also expected that P(s) ~ e~ As (A > 1) for s ^s> 1 
similar to the insulator limit P(s) = c~ s . 

A natural question to ask is whether these spectral features 
also holds at the Anderson transition in higher dimensions 
d = 5,6. Our numerical results (see Fig. 4) fully confirm 
that both P(s) and £ 2 (£) have all the signatures of critical 
statistics. The plots of P(s) and E 2 (£) correspond to the 
maximum L used in each dimension though almost identi- 
cal results are obtained for smaller volumes (not shown). The 
straight lines in Fig. 4 are fits of the form S 2 (£) — C + \l and 
In P(s) = D — As. The best fitting parameters \ an d A are 
plotted in Fig. 5 as a function of the spatial dimensionality. 
It is clearly observed that the slope of the number variance \ 
increases and A decreases with the spatial dimensionality but 
does not reach the Poisson limit x = A = 1, This confirms 
that the upper critical dimension must be d u > 6 and strongly 
suggests that it is indeed infinity as this is, according to the 
fitting in Fig. 5, the dimension in which x = A = 1. 

We are especially interested on the specific dependence of 
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FIG. 5: Dependence of the slope of the number variance \ (left 
panel) and the asymptotic decay of the level spacing distribution A 
(right panel) with the dimensionality of the space. 



X and A with the spatial dimensionality. In d = 2 + e the An- 
derson transition occurs in the weak disorder region, g c ^ 1 
and x~ l/g c ~ d—2 1. On the other hand, the prediction 
for d — * oo (Cayley tree) is A = x = I- 27 in principle, cor- 
rections to the Cayley tree limit due to interference between 
different paths decay as 1 jd or faster so it is tempting to con- 
jecture that x = 1 - C/(d - 2) and A = 1 + D/(d - 2). The 
numerical results of Fig. 5 confirm this dependence especially 
for the parameter A. In the case of x the situation is less clear. 
A reason for the discrepancy with the theoretical prediction 
could be that d ~ 3 is still far from the limit d ;§> 2 in which 
X = 1 — C/(d — 2) holds. Indeed we have observed that our 
numerical data are better described (dotted line in Fig. 5) by 
X = tanh[C (d - 2)] with C = 0.29 - 1/tt. Such a depen- 
dence of x on hyperbolic functions has already been reported 
on the generalized random matrix models 13 whose spectral 
correlations are strikingly similar to the ones at the Anderson 
transition. 

The straight lines in Fig. 5 are linear fits to the conjec- 
tured relations with fitting parameters C = 0.78 ± 0.06 and 
D = 0.55 ± 0.01. From a physical point of view these nu- 
merical results are a further confirmation that analytical ap- 
proaches to the Anderson transition starting from the metallic 
limit and adding interference corrections or starting from the 
insulator state and inducing the transition to a metal by in- 
creasing the tunneling amplitude fail to capture key features 
of the Anderson transition in intermediate dimensions where 
both mechanisms are at work. 



D. Random matrix models and the Anderson transition 

Typical signatures of critical statistics have also been found 
in both generalized random matrix model s 13 ' 32,33 whose joint 
distribution of eigenvalues can be mapped onto the Calogero- 
Sutherland model at finite temperature and phenomenolog- 
ical short-range plasma models whose joint distribution of 
eigenvalues 34 is given by the classical Dyson gas with the log- 
arithmic pairwise interaction restricted to a finite number k of 
nearest neighbors (the spectral correlations of this model are 
usually referred to as Semi-Poisson statistics though this name 
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refers to the case k — 2). In the latter explicit analytical solu- 
tions for all correlation functions are available for general k. 
Although these models reproduce typical properties of critical 
statistics such as spectral scale invariance, level repulsion and 
linear number variance, they are quantitatively different. In 
the generalized matrix models the joint distribution of eigen- 
values can be considered as an ensemble of free particles at 
finite temperature with a nontrivial statistical interaction. The 
statistical interaction resembles the Vandermonde determinant 
and the effect of a finite temperature is to suppress smoothly 
correlations of distant eigenvalues. In the case of the short 
range plasma model 34 this suppression is abrupt since only 
nearest neighbor levels interact each other. A natural question 
to ask is which of those mechanisms is dominant in the Ander- 
son transition studied in this paper. We have found a method 
to distinguish between them. In the short range plasma model 
Ax = 1 (A describes the exponential decay of P(s) ~ e~ As ). 
By contrast, in the generalized random matrix models A\ falls 
between 1/2 in the region of weak disorder to unity in the re- 
gion for strong disorder. On the other hand in our case - a 
disordered tight binding model at the Anderson transition - 
A\ ranges from 0.44 in d = 3 to 0.9 in d = 6 in agreement 
with the prediction of the generalized random matrix models. 
Our results thus suggest that the abrupt suppression of spectral 
correlations typical of Semi-Poisson statistics can describe the 
spectral correlation at the Anderson transition in d 3> 2 but 
not for intermediate dimensions. 



E. Effect of a magnetic flux 

So far all results we have presented correspond to the case 
of time-reversal invariance. We have also investigated the ef- 
fect of a random flux at criticality in d = 3 — 6. This has 
been achieved by the substitution tij — *■ tije l9ij in the Hamil- 
tonian Eq. (1). The phases 0y were chosen to be uniformly 
distributed in the interval [— tt, tt]. In d = 3, in agreement 



with previous claims in the literature- 1 , small differences with 
respect to the time reversal invariance case were found in W c 
and in P(s) in the s«l limit. Typically these effects are re- 
lated with weak localization like corrections that are strongly 
affected by the flux. However in d = 5,6 the time-reversal 
and the time-broken cases were almost indistinguishable. This 
suggests that the mechanism of localization leading to weak 
localization corrections based on destructive interference is 
less important in d 3> 2 dimensions. 



IV. CONCLUSIONS 

We have studied the dependence on the spatial dimension- 
ality of different quantities relevant in the description of the 
Anderson transition. As a result we have concluded that the 
upper critical dimension for localization is infinity. The level 
statistics tend to Poisson statistics, typical of an insulator, as 
the upper critical dimensionality is approached. We have also 
proposed that the exponential decay of the TLCF observed in 
numerical calculations is a signature of an Anderson transi- 
tion. Neither the self-consistent theory of localization exact 
in the Cayley tree nor the e-expansion formalism are accurate 
for intermediate dimensions. A new basis for the localization 
problem is thus called for. Finally, the effect of a magnetic 
flux and the validity of certain effective models to describe 
the spectral correlations at the Anderson transition have been 
discussed. 
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